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^ ■ ABSTRACT 

o ■ 

We use phase space methods to investigate closed, flat, and open Friedmann- 
D I Robertson- Walker cosmologies with a scalar potential given by the sum of 



two exponential terms. The form of the potential is motivated by the dimen- 



^ ■ sional reduction of M-theory with non-trivial four-form flux on a maximally 

H ; 

, symmetric internal space. To describe the asymptotic features of run-away 

solutions we introduce the concept of a 'quasi flxed point.' We give the 
complete classiflcation of solutions according to their late-time behavior (ac- 
celerating, decelerating, crunch) and the number of periods of accelerated 
expansion. 
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1 Introduction 

The two most intriguing (and supposedly related) recent results in observa- 
tional cosmology tell that the universe is currently undergoing an accelerated 
expansion while the main contribution to its energy density comes from dark 
energy ^1 121 El- As the origin of the dark energy remains elusive from the 
point of view of standard cosmology and particle physics p|, a number of 
possible explanations have been put forth, ranging from postulating the ex- 
istence of a non-zero cosmological constant or a scalar fleld with a suitable 
potential (a form of quintessence to modiflcations of general relativity. 
Many of these ad-hoc models are in agreement with the observational data, 
therefore it is interesting to investigate those which can be derived from a 
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fundamental, unifying theory of all interactions, a good candidate for such 
being string/M-theory. 

One particular and relatively simple approach in this direction, which has 
received considerable attention over the last year, is string or M-theory com- 
pactified on a time-dependent maximally symmetric internal manifold jE], 
optionally adding the flux of a covariantly constant four-form field strength 
[H El ini mij- Given a non-flat geometry of the internal space, the resulting 
four-dimensional low energy effective action involves a scalar field (stemming 
from the dynamical volume of the internal space) governed by a potential 
with one exponential term. Switching on flux adds another exponential term. 
This setup quite generically leads to a transient period of accelerated cos- 
mological expansion inducing a modest increase of the scale factor. In refs. 
fn\ it was pointed out that for certain values of parameters also eternally 
accelerating expansion is possible, albeit the acceleration tends to zero at late 
times. In refs. [121 CH CSl [12] these features were qualitatively inferred by 
relying on an earlier work due to Halliwell P7| based on phase space methods. 

Halliwell's original phase space analysis considered a scalar field with a 
single exponential potential in a homogeneous and isotropic universe. Subse- 
quently it was extended to include barotropic matter ^H] and multiple scalar 
fields [ini [201 [211 [22|, and generalized for anisotropic universes P3]. The 
phase space analysis of cosmologies arising from toroidal compactifications 
of string and M-theory with a four-form flux (inducing a single exponential 
potential) was carried out in refs. [^[2H]- Potentials given by the sum of two 
exponential terms have been used to model quintessence in refs. [201 [2Z1 [2H| • 
The aim of the current paper is to systematically analyze the dynamics of the 
double exponential case using phase space methods. Even though our main 
interest is in potentials arising from string or M-theory compactifications, we 
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have not limited our study to these cases and effectively cover all possible 
positive ratios of the coefficients in the exponentials. 

Phase space methods are particularly useful when the equations of motion 
are hard to solve analytically.^ Computing numerical solutions with random 
initial conditions is not a satisfying alternative as this may not reveal all 
the interesting properties. To capture the features of solutions where the 
scalar field runs away to infinity we introduce the concept of a 'quasi fixed 
point': a fixed point of the reduced dynamical system which describes the 
asymptotic limit of run-away solutions. In this manner we gain results analo- 
gous to the single exponential expected P3|, since in the asymptotic 
regime one of the exponential terms dominates over the other. New features 
occur in the regime where both exponential terms are important. These in- 
clude an unstable de Sitter fixed point which, however, cannot be obtained 
from maximally symmetric compactifications. Combining the information 
from the fixed point analysis and numerical solutions, we are able to give 
the complete classification of solutions according to their late-time behav- 
ior. In particular we discover a new class of M-theory cosmologies with two 
accelerating phases, one transient and another eternal.^ 

The paper is organized as follows. In section 2 we begin by considering 
a higher-dimensional gravitational action with a flux term and dimension- 
ally reduce it on a maximally symmetric internal space which can be either 
hyperbolic, toroidal, or spherical. We use the Friedmann-Robertson- Walker 
(FRW) ansatz with a non-trivial lapse function to obtain the equations of 

"'^The general analytic solution of flat FRW universe and a scalar field with a single 

exponential potential was recently found in ref. 29 , while for certain values of parameters 

general analytic solutions are also known when barotropic matter is added [M)\ 1811 1^ . 
^When this manuscript was in preparation, ref. appeared, mentioning a possibility 

for two periods of acceleration in this setup, but without giving specifics. 
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motion which we write in the form of a dynamical system. This helps us to 
map out the structure of the phase space (regions corresponding to accel- 
eration or deceleration and to open, flat, or closed universes) and to obtain 
detailed information about the various fixed points. With these analytic re- 
sults at hand, we discuss the phase portraits in section 3. There we also give 
a classification of all possible solutions, according to their late-time behavior 
(accelerating, decelerating, singular) and their number of phases of acceler- 
ated expansion. In section 4 we consider solutions at the fixed points and 
investigate their physical properties, e.g., the total energy density and the 
state equation parameter. We also observe that in all cases the internal space 
undergoes decelerated expansion. Finally in section 5 we conclude with some 
remarks on the viability of this setup in the light of cosmological observations 
and give a brief outlook on string theory cosmology. 



2 Analysis of the Equations of Motion 

In this section we first review the derivation of the four- dimensional effec- 
tive action obtained from the dimensional reduction of [d + m)-dimensional 
gravity coupled to m-form flux on an m-dimensional maximally symmetric 
internal manifold (see for example JH])- Taking the (i-dimensional space-time 
to be four-dimensional FRW we rewrite the resulting equations of motion as 
a dynamical system and analyze the properties of the corresponding phase 
space analytically. 
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2.1 The higher-dimensional input 

Our starting point is a gravitational theory in {d + m) dimensions coupled 
to an antisymmetric m-form field strength, 

This type of ansatz also includes M- and type II string theory with the dilaton 
consistently set to zero. We take the flux to be proportional to the volume 
form of the internal space, ^ 

Fm = bVoU. (2) 

In the case of M-theory and d = 4, Fm is a seven-form, which is dual to the 
four-form field strength which usually appears in the action. Notice, since 
this seven-form carries indices taking values in the internal space only, its 
Chern-Simons term vanishes, so that there is no obstruction dualizing to a 
four-form. 

We then make the following ansatz for the dimensional reduction to d 
dimensions: 

ds^ = ~g^,{x) dx^ dx'' + gmn{x, y) dy^ dy" . (3) 

Here g^y denotes the d- dimensional space-time metric (with coordinates x), 
and (jmn is the metric on the internal space (with coordinates y) which takes 
the form 

Qmn = e^^^""^ Cj^n{y) ■ (4) 

■^The most general flux compatible with a homogeneous and isotropic space-time has 
been considered in jl3|. where it was foimd that all other possibilities are either dual to 
this ansatz or do not contribute to the scalar potential. 
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Here e^^^^^ corresponds to a warp-factor of the internal space and cD^n is the 
normahzed metric on the maximally symmetric internal space satisfying 



Rmn{<^) = ki {m- 1) (l^mniy) , J d'^V V'^^iv) 



(5) 



The type of internal space indicated by ki can be spherical {ki — +1), toroidal 
{ki = 0), or hyperbolic {ki — —1), where the latter has been made compact 
by dividing out a suitable discrete isometry. 

Upon dimensional reduction to d = 4 and by performing a conformal 
rescaling of the space-time metric — > Qf^i,, one obtains the Einstein frame 
effective action 



(6) 



with the scalar potential 



Vi<P) ^-k, e-2^<^ + ^e-?^ . (7) 

Here we have set the gravitational couplings to unity, k^^^ — 1, = 1, 

and introduced c = \J^^^- The quantities h and (j) ^-^e related to their 
microscopic counterparts h and [3 by 

62^:46M^ ] , (8) 

Vm m- 1)7 ' ^ ' 



= ct>{x) + ^— In ^ , (9) 



and 

fl(x^ = , 

mc m + 2 \ 4 

respectively. 

By dimensional reduction we can obtain 1 < c < a/3 only, where the case 
c = -\/3 corresponds to compactification on S^. Nevertheless we will analyze 
the full interval < c < oo for two reasons, first for the sake of completeness 
and second because other types of compactifications might lead to c-values 
outside the interval 1 < c < \/3. 
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2.2 Four-dimensional FRW cosmology 

Let us now consider FRW space-times with a non-trivial lapse function, 

ds^ = -N{tYdt^ + e^-W ( ^-—dr^ + . (10) 

\1 — kr^ J 

Here denotes the line element of the two-sphere and k = —1, 0, 1 corre- 
sponds to an open, flat, or closed universe, respectively. Using this ansatz 
we obtain Friedmann's equation, 

as well as the dynamical equations, 

where the dot denotes a derivative with respect to the time variable t. The 
cosmological time r is related to t via dr = N{t)dt. 

In order to study these equations as a dynamical system we proceed as 
follows. First we use Friedmann's equation to eliminate the term containing 
a{t) from the dynamical equations. Then we specify the lapse function N{t): 

N{t) := ^e'<^W, fci = 0, 

N{t) := e^<^W, ki = ±l, c < , (14) 

N{t) := ^ei<^W, h = ±1 , Vs<c. 

For finite values of (p{t) these lapse functions are positive definite, implying 
this is an admissible choice. As long as no approximations are made different 
lapse functions lead to different parametrizations of the same physics. The 
choice ()14|1 is motivated by the regime (p oo, where one of the exponentials 
appearing in the potential ((Zj) becomes unimportant and may be truncated. 



In this approximation the effect of the leading exponential term is encoded 
in the definition of our time variable. As explained below eq. ()40|) this 
procedure gives the correct late-time behavior of the solutions. The case 
c = y/S, fci = ±1 is special, as for this particular value of c the exponentials 
appearing in the potential are the same, V = — ki^ g-2\/3(/> Therefore 
this case can be treated analogously to the single exponential case, ki = 0, 
and will not be considered explicitly. 
Finally we introduce the variables 

x{t) := (Pit) , y{t) := , z{t) := a{t) , (15) 

in order to rewrite the dynamical equations (O, (fT^ as a set of autonomous 
first order differential equations. In the single exponential case {ki = 0) we 
obtain the two-dimensional dynamical system I^'(y, z): 

{ y= 3_ 2 _^^^G 

Note that we do not take the trivial equation a; = ?y to be part of the dy- 
namical system. This is possible because the equations for y and z do not 
depend on x, as a result of the choice of the lapse function (fT^ . Therefore 
we have an autonomous system which only involves y and z. This system 
has been studied in refs. jT3 IT^ ITT)]. 

In the double exponential case {ki = ±1) we obtain the three-dimensional 
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dynamical system V{x,y, z) 



X = 








y = 


cy2 - 


-3yz- 




z = 


cyz 




z'-lh + '^ 


X = 


y, 






y = 




-3yz + 


C fe2 


z = 


lyz 




2,2 iki 
^ ^ 3 3P 



3 e 



^-3)x 
-3)x 



< Vs. 



(17) 



In this case we cannot drop the equation x = y, because the equations for 
y and z depend on x. The x-dependent terms correspond to the subleading 
exponential term, while the leading one, which dominates for large (j), was 
absorbed by the choice ([1^1 of the lapse function. 

2.3 The structure of the phase space 

After defining the variables parametrizing the phase spaces of our systems, 
we will now discuss their structure. As the dimension of the phase space for 
A;i = and fci = ±1 is different, we treat them separately in the following. 

We start by investigating which points in phase space are associated with 
curvature singularities. To this end we compute the Ricci-scalar of the metric 

R = ^(^a + 2a^-^a + kN^e-^'^^ . (18) 

Substituting in Friedmann's equation, the specific choice of the lapse function 
()14|). and the dynamical equation ()12|) this expression takes the following 
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form: 

R = 62g-fx (4_2;2) ^ h = 0, 

R = 2e-2- (^-4fci + 2 62ei(^'-=^)"-y2j ^ h = ±1 , c < Vs , 

(19) 

These equations indicate that for both ki = and ki = ±1 the boundaries of 
the phase space x = — oo and y = ±oo correspond to curvature singularities, 
while at the surfaces z = ±oo as well as x = +oo the curvature scalar remains 
finite. This result also holds for other curvature invariants as i?^,^ R^'^ and 

Next we determine the phase space regions, where the solutions undergo 
acceleration in the scale factor a(T) = e"'-*'-'^^'' . These are characterized by 
the inequality 

„M":.i^„M.(^^)%«(..>0, (20) 

Substituting in the lapse functions (fT^ together with the dynamical equa- 
tion (fT^ . this implies that acceleration occurs for 

2/2 < 1, A;i = 0, 

y2 < _fc^ + |el(c^-3)-^ fci = ±l, c<v^, (21) 
2/2 < i_^e-K-'-3)x^ A;i = ±l, V3<c. 

Likewise the regions of expansion are given by ^ |e"W > 0, which corre- 
sponds to 2; > 0. 

From eq. ^2\\ we observe that acceleration occurs if 0^ (and hence the 
kinetic energy of the scalar field) is small. We further find that the flux 
contribution enlarges the acceleration region if c < -\/3 or /ci = — 1 while for 
-\/3 < c, fci = +1 the flux narrows the acceleration region. 
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Finally, we can use Friedmann's equation to separate the phase space into 
regions corresponding to closed, flat, or open universes. Setting A; = and 
substituting in eq. (|T^ we obtain the k = hypersurfaces: 

-3z^ + 2 = 0, ki = 0, 

y2_Q^2_2k, + pe-^i'^'-^)- = 0^ k^ = ±l,c<VS, (22) 
2/2 - 3 ^2 + 2 - Ml e-K-'-3)- = , ki = ±l,V3<c. 

For the left hand side being smaller (bigger) then zero, we are in the k = —1 
{k = +1) region, respectively. 

2.4 Fixed points 

Let us now determine the fixed points (critical points, equilibrium points) 
of the dynamical systems T>'{y, z) (fTBj) and V{x, y, z) (fTTj) . Recall that for a 
single exponential potential {ki = 0) it is possible to eliminate both a (using 
Friedmann's equation) and x = </) (using the lapse function ()14|l ). so that the 
dynamical system is two-dimensional. The fixed points (FPs) are given by 

y Ipp = , i Ipp = . (23) 

At a FP the scalar field (j){t) and the scale factor a{t) therefore have the 
following form: 

0(t) = 0* t + ci , a{t) =a*t + C2. (24) 

Note that <p does not go to a constant at the FPs of V{y,z). That would 
be the case if we looked for FPs of the larger system V{x,y,z), which is 
obtained by adding the third equation x = y to (fTBj) . However, for a single 
exponential potential V{x,y,z) does not have FPs. This is clear, because 
the potential does not have an extremum. The generic late-time behavior of 
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eternal solutions is a run-away behavior ^ oo with asymptotically constant 
velocity (with respect to the time variable t), corresponding to the FPs of 
V'{y,z). 

For of the double exponential potential it is not possible to eliminate 
X = (p from the equations for y and z and we have to deal with the three- 
dimensional dynamical system V{x,y, z) (fTTj) . The FPs of this system satisfy 
the equations 

2/ Ifp = , i; Ipp = and x |pp = , (25) 
and solutions for 0(t) and a{t) at the FP take the form 

(P{t) = (p*, a{t) = a*t + ci. (26) 

A double exponential can in principle have an extremum at a finite value of 0, 
and we will see later that for a certain range of parameters the corresponding 
FP indeed exists. 

We also expect that solutions with the run-away behavior (j2l|) exist and 
play a role in the dynamics of V{x,y,z). Observe that for — oo the x- 
dependent terms in ()17|) become exponentially small. Therefore it is a reason- 
able approximation to drop them, which leads to a reduced two-dimensional 
dynamical system T'oo(?/, z).^ For ^/3 < c this is in fact the same system as 
the one obtained for ki = 0, T>'{y,z) (fTI)|) . The physical picture behind this 
truncation is that solutions which move towards (p = oo enter an intermedi- 
ate regime of large but finite (p, where the subleading exponential becomes 
irrelevant and can be dropped, while the leading exponential is encoded in 
the lapse function (fT^ . The behavior in this regime is then controlled by 

the reduced dynamical system Voo{y,z), i.e., cosmological solutions adjust 
'^If one just drops the first equation in (|17() the resuhing non-autonomous system does 
not have FPs for finite x. 
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to the FPs of Voo{y, z) (or are repelled from, depending on the nature of the 
FP). In the three-dimensional parameter space {x,y,z) we expect that the 
solutions flow in a narrow tube towards the attractive FPs of Voo{y, z), which 
are located at the boundary x = oo. The FPs of Voo{y,z) will be referred 
to as quasi fixed points (QFPs) of the full system V{x, y, z) in the following. 
At such QFPs, the time-dependence of (f){t) and a{t) is the same as the one 
found for fci = 0, ()24|1 . For terminological convenience we will also refer to 
the FPs of the reduced system z) as QFPs.^ 

We remark that the above definition of a QFP is tied to a particular, 
physically motivated approximation, namely dropping the sub-leading expo- 
nential. The effect of the leading exponential is captured by the lapse func- 
tion, which therefore acquires a physical meaning. Comparing to numerical 
solutions of the full equations of motion in the next section, we establish that 
the picture outlined above is indeed correct and that the concept of QFPs is 
useful for describing the dynamics of y, z). Also note that the behavior 
of (j){t) and a{t) at a QFP of y, z) involves a distinguished time coordi- 
nate. In contrast to a true FP in phase space, ()24|1 is not reparametrization 
invariant. Of course, reparametrization invariance as such is not an issue, 
since we use a specific (class of) parametrization(s), ()10j). anyway. However, 
the time variable t does not have a direct physical interpretation. Therefore 
we will investigate the analytical properties of (Q)FP solutions with respect 
to the cosmological time r in section 4. 

Whether a FP or QFP is an attractor, in the sense that there is a class of 
solutions which asymptotically evolve with the same values (and (^*), 

can be deduced form the eigenvalues of the stability matrix 

Bijipp dj(3i Ipp . (27) 

^Note that the FPs discussed in [13 ESI ESI are QFPs in our sense of definition. 
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Here jSi is short for the right hand sides of eqs. p6|) and (fTTj) and i G {y, z}, 
X = [y, z]^ , or z G {x, y, z}, X = [x, y, z]^ , depending on whether we evaluate 
the condition ()23j) or ()25|) . Denoting the eigenvectors and eigenvalues of this 
matrix as Va and 6a, respectively, we can write down the linearized solution 
of the full system in the vicinity of a (hyperbolic) FP, 

X{t) = X* + J2caVae'-\ (28) 

a 

with Cq being the constants of integration. This indicates an attractor has 
9a < for all eigenvalues, while for mixed positive and negative eigenvalues, 
one obtains a saddle point. 

An interesting subtlety in the attractor behavior arises if the (Q)FP is 
a saddle point located on the k = surface. In this case it is necessary to 
check the direction of the eigenvectors relative to the normal of the surface 
to deduce for which solutions the (Q)FP acts as an attractor. For z > 
(expansion region) the typical configuration arising here is that the attractive 
eigendirections are tangent to the surface while the repulsive ones have a 
non-zero component along the normal vector, such that the (Q)FP acts as 
an attractor for A; = solutions while solutions starting in the regions k = ±1 
are repelled and cannot reach it. For z < the opposite situation occurs. 
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The case ki = 

We start by analyzing the QFPs of eq. pfij) . Exphcitly evaluating the QFP 
condition pHjl yields that the system has two pairs of QFPs:^ 

QFPi = [±1 , ±|] , exists for < c < oo , 

(29) 



QFP 



2 



c2-3 ' 3 V c2-3 



exists for -\/3 < c . 



Substituting the QFP coordinates into the equations for the acceleration re- 
gion (PT|) and the k = surface yields that QFP^ is located in the 
k = —1 and A; = +1 region for < c < 3 and 3 < c, respectively. Further it 
is always situated at the boundary between the acceleration and deceleration 
region. The QFPg always resides on the A; = surface, being in the deceler- 
ation region for V3 < c < 3 and in the acceleration region for 3 < c. The 
case c = 3 is special. Here the two QFPs merge into a single non-hyperbolic 
QFP at [±1,±1]. 

To determine the properties of the phase space trajectories close to the 
QFPs, we calculate the eigenvectors and eigenvalues of the stability matrix 
(j2Z|): 



= i (-3 + V4c2 - 27) , = i (-3 - V4c2 - 27) , 



^2+ = h/^(6-i^^)> ^yS(3-c^)• 



(30) 



These eigenvalues imply that the QFP^ is a complex attractor {i.e., the 
solutions spiral into this point), a real attractor, and a saddle point for < 
c < |-\/3, |-\/3 < c < 3 and 3 < c, respectively. Analogously one obtains 
that the QFPg is a saddle point for ^/S < c < 3 and becomes an attractor 
^Due to the invariance of (|16|l under time-reversal, (Q)FPs always come in pairs, i.e., if 
[(p* , a*] is a (Q)FP, so is [—4>*, —a*]. Further the eigenvalues of the second ones are given 
by the negative of the eigenvalues of the first, so that in the following we mention the first 
ones only. 
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for 3 < c. These results are also summarized in tabled Note that this table 
includes the QFPs in the expansion region [z > 0) only. 

In the regime ^/3 < c < 3 the QFP2 is a saddle point on the k = 
surface. By comparing the eigendirections of the QFP with the normal of 
the surface we find that the attractive eigendirection is tangent to the A; = 
surface, while the repulsive eigendirection has a non-trivial component along 
the normal. This indicates that the QFPg acts as an attractor for solutions on 
the k = surface while trajectories starting in the k = ±1 regions are repelled 
and cannot approach it asymptotically. These features are also illustrated in 
the second diagram of figure [H 

The case ki = ±1 

Let us now consider the case fci = ±1. We start by analyzing the condition 
for 'true' FPs. Imposing the conditions ()25|1 . the only non-trivial solution 
requires fci = +1 and is given by 

exists for c > . 

(31) 

It is located on the k = surface and inside the acceleration region. In 
fact, it will turn out in section 4 that this FP gives rise to an exponentially 
growing scale factor in cosmological time and therefore is a true de Sitter FP. 
Taking i G {x, y, z} and calculating the eigenvalues of the stability matrix 
(P7|l for this FP yields 

91 = V2(c2 - 3) , 9l = -2^2{c^-3), Ol = -'^ ^2{c^ - 3) . (32) 

Hence it is a saddle point having two late-time attractive and one late-time 
repulsive eigendirection. Comparing its eigenvectors to the normal of the 
k = Q surface yields that one attractive and the repulsive eigendirection are 
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(Q)FP 


h 


c 


type 


k 


acceleration 


1 





< c < 1 V3 


complex attractor 


-1 


boundary 






1 V3 < c< 3 


real attractor 


-1 


boundary 






3 < c 


saddle point 


+1 


boundary 


2 





< c < V3 


— 


— 


— 






V3 < c < 3 


saddle point 





deceleration 






3 < c 


real attractor 





acceleration 


3 


+1 


< c < V3 












< c 


saddle point 





acceleration 


4 


-1 


< c < 1 


saddle point 


+1 


boundary 








real attractor 


-1 


boundary 






^<c<V3 


complex attractor 


-1 


boundary 






V3 < c 


— 


— 


— 


5 


-1 


< c < 1 


real attractor 





acceleration 






1 < c < V3 


saddle point 





deceleration 






< c 


— 


— 


— 


6 


±1 


< c < 


— 


— 


— 






V3 < c < |V3 


complex attractor 


-1 


boundary 






|V3 < c < 3 


real attractor 


-1 


boundary 






3 < c 


saddle point 


+1 


boundary 


7 


±1 


< c < V3 












< c < 3 


saddle point 





deceleration 






3 < c 


real attractor 





acceleration 



Table 1: Summary of the fixed point properties fomid in subsection 2.4. Here '(Q)FP' 
gives the number of the (quasi) fixed point as it appears in the text, while and 'c' are 
the parameters of the potential The columns 'type', 'A:' and 'acceleration' specify the 
properties of the (Q)FP for the corresponding values of the parameters. At the values of 
c not covered in this table we encounter non-hyperbolic FPs. 
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tangent to the surface while the second attractive eigenvector has a non- 
zero component along the normal. This implies that the FP3 is a saddle 
point on the k = surface, i.e., flat cosmological solutions need to be fine- 
tuned in order to approach it asymptotically. Additionally, we can also fine- 
tune solutions starting in the k = ±1 region such that they asymptotically 
approach this FP. The latter properties are reflected in table El while the 
results on the position of the FP are compiled in table ^ 

Notice that the FP3 is located at the (positive) maximum of the potential 
(|7|) with ki = +1. This explains why it is unstable (saddle) and corresponds 
to an expanding de Sitter universe. For c < a/3 the potential has a (negative) 
minimum, but there is no related FP, since this situation corresponds to 
an anti-de Sitter universe which contracts to a crunch represented by the 
boundary of the phase space. For ki = —1 the potential has no extrema and 
we do not obtain any FPs in this case. 

Having found the true FP of the system (fT7|) we now relax the FP con- 
dition and look for QFPs satisfying (j2SI)- As previously observed, the QFPs 
are always located at the x = 00 boundary, but fix the values <f)*,a* along 
which the boundary is approached. 

In the case c < -\/3, where the ki term is dominant, substituting the 
ansatz (j^^ into (fTTjl and dropping the subdominant exponential term only 
gives rise to QFPs for ki = —1: 




QFP4 = [±1 , ±c] 



exists for c < -\/3 , 
exists for c < \/3 . 



(33) 



The corresponding eigenvalues of the stability matrix are given by 



e^_^ = -c + V4 - 3 c2 , 



-c - V4 - 3 c2 



4 




-V2(3-c2). 



(34) 



'5+ = 
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The properties of these QFPs with respect to acceleration and the A; = 
surface are, again, obtained by comparing the position of the QFP relative 
to the surfaces ^IT\ and (j221), see table [U 

In the case c > y/S the flux-term gives the dominant contribution to the 
potential and substituting x = oo into (fTTj) is equivalent to setting ki = 0. As 
a consequence, we obtain the QFPs pHj) together with the properties (jSDI)- 
For completeness these QFPs are labeled QFPg and QFP^, respectively, and 
their properties are also included in table H 

3 The Phase Portraits 

The greatest advantage of the phase space method is that one can plot phase 
portraits, which makes the classification of solutions easy and intuitive. Here 
one uses two well-known properties of first order ordinary differential equa- 
tions: (i) trajectories (integral lines) do not cross, and (ii) for hyperbolic 
FPs the local behavior of the exact solution is uniquely determined by the 
linearized system. The first point can often be used to show that the phase 
space decomposes into disjoint subspaces. The second point helps the analy- 
sis, because the hyperbolic QFPs found in the last section govern the behavior 
of the trajectories in their vicinity and lead to an universal asymptotic be- 
havior. We will first explain the two-dimensional phase spaces occurring in 
toroidal compactifications with flux. These results are already covered by 
the literature ITK | ITI)] . but serve well for explaining the method. Then we 
turn to the more complicated phase spaces of curved flux compactifications. 
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3.1 Toroidal internal space with flux 

The properties of the two-dimensional phase spaces arising in models with 
ki = are summarized in figure Q and table El They include all fiux com- 
pactifications with toroidal internal space. The QFPs of this system were 
classified in the last section and their properties are summarized in tabled 
Based on this classification, there are four different cases for the parameter 
c. Note that the classification of hyperbolic compactifications without fiux 
is also covered by table El and requires a redefinition c — > |. 

Figure Q shows a representative phase portrait for each of the four cases. 
Here the horizontal and vertical axis correspond to y = (p and z = a, respec- 
tively. The universe expands for z > and contracts for 2; < 0. These two 
regions are mirror images of one another as they are related by time reversal, 
which acts on V{x,y,z) hj x —>■ x, y ^ —y and z —z. The vertical 
gray strip specifies the region of acceleration determined by eq. ()21|) . More 
precisely, this region corresponds to accelerated expansion for z > and de- 
celerated contraction for 2; < 0. Eq. implies that, independent of the 
value of c, the trajectories corresponding to fiat universes (black lines) cut the 
phase space into three horizontal regions, corresponding to open universes 
(expanding in the upper and contracting in the lower region) and closed 
universes (middle region). The specific late-time behavior of a trajectory is 
encoded in its color and can be inferred from table El 

As an example, let us consider the phase portrait in the upper left corner 
of figured It corresponds to the M-theory value c = In this range of 
c only QFPi exists, coming in two copies: complex attractor for 2; > and 
complex repeller for z < 0. These are located in the regions of open uni- 
verses and at the boundary between acceleration and deceleration. All open 
universes in the upper {z > 0) region start at a boundary and spiral into the 
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Figure 1: Phase portraits for fci = 0. The flow is shown for the four representative 
cases c = (upper left), c = 2 (upper right), c — 2.8 (fower left) and c = 4 (lower 
right), respectively. The thick black lines represent the fc = surface and the trajectory 
corresponding to a flat universe. The arrows point into the direction of increasing time 
variable t. 
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Table 2: Colors used in figures H and El to distinguish cosmological solutions according 
to their number of phases of accelerated expansion and late-time behavior. 



QFP for late times (cyan trajectories). Correspondingly, in the lower {z > 0) 
region the open universes start at the QFP and spiral out to the bound- 
ary (magenta trajectories). Flat and closed universes run from boundary to 
boundary thereby crossing the acceleration region once. Since the expan- 
sion and contraction parts of the phase space for closed universes are not 
separated, the behavior of solutions can change along the trajectory: in the 
acceleration region contraction can become expansion while in the decelera- 
tion region expansion can switch into contraction, but not vice versa. Hence 
the closed universes can have either one (violet) or no (magenta) phases of 
accelerated expansion. In summary we get two types of late-time behavior 
possible in the whole setup: expanding open universes exhibit eternal oscilla- 
tions between acceleration and deceleration, while contracting open universes 
as well as all flat and closed universes terminate at a boundary of the phase 

space associated with a curvature singularity.^ 

^Numerical investigations of the behavior of the logarithmic scale factor for closed 
universes indicate that the boundary is reached in finite cosmological time and corresponds 
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Table 3: Classification of possible cosmological solutions for fci = according to their 
late-time behavior and their number of periods of accelerated expansion. Here 'G' denotes 
that the corresponding behavior is generic, 'F' indicates that the behavior requires a fine- 
tuning of the initial conditions, and 'C marks a special solution which runs from one fixed 
point into another. 



Similar results can be read off for the other ranges of c displayed in figure 
The main addition is that we do also encounter examples of special, non- 
generic solutions. First, there are fine-tuned solutions which run into or start 
from a saddle point and are thus future or past infinite, respectively. Second, 
there are so-called cross-over solutions which run from one fixed point into 

another and thus have both infinite past and infinite future, 
to a —> — oo, i.e., a big crunch. 
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The classification of all possible solutions is given in table El which is 
organized as follows. The first and second column give the values of c and k. 
The following columns characterize the late-time behavior: singular, meaning 
that the solution runs to the boundary of the phase space, and accelerated 
expansion (a" > 0, a' > 0), decelerated expansion (a" < 0, a' > 0), acceler- 
ated contraction {a" < 0, a' < 0), decelerated contraction {a" > 0, a' < 0), 
or cyclic (oscillation between accelerated and decelerated expansion) for so- 
lutions asymptotically approaching a QFP. Since some solutions also exhibit 
transient accelerated expansion, we have split the columns further to indicate 
the total number of phases of accelerated expansion. (For the cyclic case this 
number is infinite). Within the table a dash indicates that the correspond- 
ing behavior does not occur, while 'G' says that it is generic. By generic we 
mean that one does not need to tune initial values to find a trajectory in 
this class. In contrast, 'F' represents fine-tuned solutions, where one starts 
from a locus of codimension at least one in the space of initial conditions. 
Finally 'C stands for a solution which "crosses over" from one fixed point to 
another. 

3.2 Curved internal space with flux 

We now turn to the three-dimensional phase spaces of models with ki = ±1, 
which include, for suitable values of c, flux compactifications on curved inter- 
nal spaces. Although the method remains the same, one faces new practical 
challenges here. Three-dimensional phase spaces are harder to visualize, es- 
pecially in printed media which only allow static two-dimensional projections. 
Therefore we have only included one "three-dimensional" plot, figure El for 
illustrational purposes, and summarize our results in tables E] and El 

For the classification of the trajectories according to their accelerational 
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Figure 2: Phase portrait for hyperbolic flux compactifications of M-theory, c — ^ (and 
taking 6=1). The diagram shows one example for each class of trajectories encountered in 
the k — —1 region. From right to left: eternal deceleration, eternal deceleration preceeded 
by a transient accelerating phase, eternal acceleration and eternal acceleration with an 
earlier transient accelerating phase. 

phases, it turns out that projecting the full solutions onto the (x, y)-plane is 
a very useful tool. This is due to the equations describing the acceleration 
region (pT|l being independent of z, so that the number of accelerating phases 
can directly be read off from this projection. In order to determine which 
regions of initial conditions may give rise to different behaviors of solutions 
it is further helpful to study the vector field generating the trajectories. To 
find the fine-tuned solutions, knowledge about the (Q)FP structure of the 
phase space is crucial. 

Let us illustrate the properties of cosmologies with a three-dimensional 
phase space using the most interesting examples, expanding open universes in 
hyperbolic flux compactifications of M-theory, see figure |21 Here the horizon- 
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tal and vertical axis correspond to y = <p and z = a as before, while the new 
direction orthogonal to this surface is the new phase space variable x = <p. It 
is useful to compare with the phase portraits of figure ^ which correspond 
to the X ^ +00 limit of the three-dimensional plot. As is apparent from this 
comparison, the two surfaces aligned to the (x, y)-plane are the subspaces 
corresponding to flat universes. Thus again we have disjoint regions popu- 
lated by expanding open, closed, and contracting open universes. The two 
vertical surfaces enclose the region of acceleration between them. This is the 
generalization of the gray strip of figure H According to table ^ two (pairs 
of) QFPs exist for fci = — 1, 1 < c < namely QFP4 and QFP5. These are 
actually located at x = = 00, but we have represented their projections 
into the displayed region of the phase space by little crosses. The QFP4 is a 
real attractor in the open universe region, located on the boundary between 
acceleration and deceleration, while the QFP5 is a saddle point on the surface 
of flat universes, located in the deceleration region. It is attractive for fiat 
universes, but repulsive normal to the k = surface. Therefore expanding 
open universes approach the QFP4 with eternal deceleration or acceleration, 
depending on from which side they come. 

In order to analyze whether periods of transient accelerated expansion 
exist, one must plot and analyze the phase space flow. We found four different 
types of generic behavior: eternal deceleration (dark brown line), eternal 
deceleration with one period of transient acceleration (green line), eternal 
acceleration (blue line), eternal acceleration with one period of transient 
acceleration (red line). The existence of generic M-theory solutions with two 
periods of accelerated expansion (red line) is a new feature which does not 
exist for a single exponential potential. Therefore one needs both flux and 
a curved internal space to have such behavior. We also note that figure El 
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clearly illustrates that the concept of a QFP introduced in the last section is 
indeed useful as the trajectories flow to a; oo at the values [(f)* , a*] of the 
QFP4. 

The classification of the possible trajectories in the other c regimes and 
k regions proceeds completely analogously. Our results for hyperbolic {ki = 
— 1) and spherical {ki = +1) internal space with flux are summarized in tables 
m and El respectively. An interesting feature with spherical compactifications 
is that the k = surface becomes connected, so that the regions of expanding 
and contracting open universes are not separated any more. In this case there 
are fiat and open universe solutions which start in the 2; > region and have 
a period of transient accelerated expansion, and then go down to the 2; < 
region, where they display a period of transient decelerated contraction.* 
The FP3 occurring for ki = +1 and c > Vs has a special status as solutions 
that evolve at the FP correspond to true de Sitter solutions. This behavior is 
denoted by the subscript 'ds' appearing in table |S1 It is only available upon 
fine-tuning, as the solution has to start at this FP. The two copies of FP3, 
one in the 2: < region and one in the z > region, are connected by a 
special cross-over solution which runs through the k = +1 domain from the 
lower to the upper FP3. 

Finally let us make a remark about this classification. One should be 
aware of the fact that our results depend on a mixture of analytical and 
numerical results. Knowledge about the (quasi) fixed points basically deter- 
mines what kind of late-time behavior is possible - the solutions either run 

to a (Q)FP or the boundary of the phase space. Numerical investigation 
*Since our models satisfy the null energy condition, flat and open cosmologies cannot 
switch from contraction to expansion, but the opposite is allowed. Closed cosmologies are 
not constrained. 
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Table 4: Classification of possible cosmological solutions for fci = — 1 according to their 
late-time behavior and their number of periods of accelerated expansion. 
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Table 5: Classification of possible cosmological solutions for fci = +1 according to their 
late-time behavior and their number of periods of accelerated expansion. The subscript 

'ds' indicates the existence of true de Sitter solutions if we fine-tune our initial conditions 
to be at the FP3. In addition there are also fine-tuned and cross-over solutions which 
approach this fixed point at late times. 
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is needed mainly to determine how many periods of accelerated expansion 
can occur. Consistency checks^ indicate that the numerical results are trust- 
worthy if we evolve trajectories with initial conditions in the 2; > region 
forward in time. However, when following trajectories back in time we find 
indications that the numerics becomes untrustworthy below a certain criti- 
cal time tc where the numerical value of tc depends on the initial conditions 
chosen. Due to this fact it is hard to make definite statements about the 
early-time behavior of our trajectories. 

4 Physics at the Fixed Points 

In this section we use the information about the (Q)FPs obtained in section 
2.4 to investigate the properties of solutions asymptotically approaching these 
points. 

Fixed point solutions in cosmological time 

In order to find the physical interpretation of our solutions we need to go back 
from the time coordinate t, which was introduced to simplify the analysis of 
the equations of motion, to the cosmological time r. Our starting point is eq. 
((2H1)- At a (Q)FP, this equation is easily integrated and yields the following 
solutions parametrized by the lapse function time t: 



^One example of an internal consistency check is to compute fc as a function of the 
fields. Analytically we know that k takes the discrete values 0, 1, —1, and when plugging 
in numerical solutions deviations should be small. 



a{t) = a*t + C 
a{t) = a*t , 



(f){t) = (f)*t, for the QFPs , 
(f){t) = (j)* , for FP3 . 



(35) 
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This is the most general solution of eq. (j28j) as in the first case the constant 
of integration appearing in the (p{t) equation can always be absorbed by a 
shift of the time variable t. In the case of FP3 we have used this freedom 
to absorb the constant appearing in the a{t) equation, since (p* cannot be 
set to zero. Concerning the remaining constant of integration, (, we have to 
distinguish between the cases where the QFP is located on the k = surface 
or in the k = ±1 region. In the former case ( is indeed an undetermined 
constant, while in the latter case it is fixed by Friedmann's equation: 

e'' = ~Jl'\, forQFPi,6, = for QFP4 . (36) 

(c2 - 9) 1 - 

Switching to the cosmological time r via dr = N{t)dt with N{t) being 
the lapse function given in eq. p4|l . we can use the explicit expression for 
(p{t) at the (Q)FP to calculate the relation between the lapse function time 
t and the cosmological time r: 

^ = d^l^^-) forQFP,,,e,v, 

t = ^ In (^c0*r) for QFP4,5 , (37) 

t = j^r for FP3 . 

Using these relations it is now straightforward to write down the cosmologi- 
cal solutions at the (Q)FPs parametrized in cosmological time. Going back 
to the scale factor air) = e"*-**-'^-''' we find 



QFPi,2,6,7 : <t) = (^) e< , m = ^ In r) , (38) 

QFP^^s : a(r) = [c 0*) , 0(r) = - In (c 0* r) , (39) 



FP3 : a(r) = exp rj , 0(r) = 0* . (40) 
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Note that the evolution of the metric and the scalar fields as functions of 
r is well behaved and universal for r > 0. This confirms that our choice 
of lapse function and the notion of a QFP (which depends on the choice of 
lapse function) is physically meaningful. Switching between t and r does 
not affect the qualitative late-time features of the solutions: they are either 
governed by the (Q)FP which in cosmological time corresponds to the asymp- 
totic behavior given above, or end up in a big crunch singularity which is 
a coordinate-independent notion. This assures that the classification of so- 
lutions with respect to t carries over to the physically relevant cosmological 
time r. 

Properties of the fixed point solutions 

All QFP solutions show a run-away behavior of the scalar field, (p ~ Inr, 
which slows down to a stop as r ^ oo. The scale factor follows a power 
law expansion a ~ r'^, where the "critical exponents" A depend on the 
model parameter c only. The values of A at the (Q)FPs are summarized 
in the column 'a' of table IHl Accelerating expansion corresponds to A > 1 
and hence requires c > 3 for QFP2,7 or c < 1 for QFP5. These, however, 
are outside the range of compactifications. From compactifications we may 
obtain at most A = 1. This corresponds to a QFP on the boundary of the 
region of acceleration and deceleration. In the previous section we found that 
such a QFP can be approached from the acceleration region. In this case the 
acceleration is eternal but approaches zero so fast that one does not obtain 
a significant growth of the scale factor. This was also observed in [TTl IT^ . 

Besides the QFPs we have also the de Sitter fixed point FP3, which does 
not exist for potentials with one exponential term. Here the scale factor 
grows exponentially, a ~ expr, while the scalar field has a finite value (p*. 
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However, this FP does not exist for compactifications: (f)* goes to zero for 
c — > cxD and moves off to infinity when approaching the maximal value for 
compactifications c = ^/3 from above. Moreover, it is a saddle point. Hence 
eternal expansion requires an infinite fine-tuning, which means to pick exact 
initial values. Generating any significant growth of the scale factor by passing 
near the saddle point still requires a drastic fine-tuning. 

We can also compute the density parameter Q and the equation of state 
parameter w for the (Q)FP solutions. Recall the standard expressions for 
the density and pressure of a scalar field, 

p=/+2l^(0), p=/-2VicP), (41) 

where the prime denotes the derivative with respect to cosmological time 
T. Defining the Hubble parameter H := ^ and introducing the density 
parameter Q := Friedmann's equation (|TT|) can be written as 

Q = l + ^. (42) 
a 

Here = 1 is the critical density corresponding to a flat universe, while 1] > 1 
and Q < 1 correspond to a closed and open universe, respectively. Taking 
the k values for the (Q)FPs given in table Q together with the solutions for 
a(r) obtained by evaluating the eqs. (jHHj) . (jH^ and pn|) . we can use this 
relation to find the density parameter Q at the (Q)FP. The results depend 
on the parameter c only and are summarized in the column denoted 'fi' of 
table El 

Furthermore, the density and pressure of the scalar field are related by 

an equation of state, p = wp, where p and p are given in (jUJ. At a (Q)FP 

w takes a constant value, which is an intrinsic property of the (Q)FP, in the 
-'^"A numerical solution placed directly on FP3 with 10 digits precision triggers 10 e-folds 
of expansion before leaving the FP due to accumulated numerical errors. 
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sense that it depends on the model parameter c only. The values are again 
listed in table IHl In the compactification range 1 < c < v^, we typically 
obtain A < 1, w > -| and ^] = C(l). There is a close relation between 
the stability properties of the (Q)FP and its w value. While w < — | always 
corresponds to an attractor, w > characterizes a saddle point. For the 
boundary value w = — | the (Q)FP is attractive when located in the k = —1 
region (open universe) while it is a saddle point if k = +1 (closed universe). 
An Q = 1 attractor sitting at the k = surface requires w < A > 1 and 
is only realized for fci = — 1, c < 1 and c > 3. 

Decompactification of the extra dimensions 

Except at the saddle point FP3, the scalar field cannot be stabilized and 
exhibits run-away behavior. For the values 1 < c < -\/3, which correspond 
to compactifications of higher-dimensional theories, this implies the decom- 
pactification of the internal manifold. In order to see whether such behavior 
is in immediate conflict with observations, we need to characterize this quan- 
titatively. From eqs. Q - ® we obtain 

Mnt = j d^y = e™'^^*) j (TyV^ = e""^^'^ , (43) 

which shows that the radius of the internal manifold behaves as Tint ~ e^*^*-*. 
Taking the definition of (3{t) given in eq. and substituting the solutions 
()38|1 and (j39|l we find a polynomial increase of Tint ~ t"^, where q < 1 for 
all cases. The explicit values of q are given in the column Vint' of table IHl 
The result g < 1 implies that the expansion of the internal space always 
decelerates. 

M-theory example 

As an explicit example let us consider the hyperbolic flux compactification 
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Table 6: Summary of the properties of cosmological solutions parametrized in cosmolog- 
ical time r at the (quasi) fixed points found in section 2.4. 
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of M-theory (m = 7, c = From table IHl we find that there exist solutions 
attracted to QFP4. At this QFP the universe expands at a constant rate 
a ~ T, while the radius of the internal manifold exhibits decelerated growth, 

2 

'"int ~ ■ The energy density and equation of state parameter have the QFP 
values = ^ ~ .77 and w = respectively. 

5 Discussion and Conclusions 

In this paper we used phase space methods to investigate the properties of 
Friedmann-Robertson- Walker cosmological solutions with a scalar field gov- 
erned by a potential with two exponential terms. The form of the potential 
is motivated by the compactification of M-theory on a maximally symmetric 
internal manifold including flux. Our main result is the complete classifi- 
cation of all possible open, flat, and closed cosmologies according to their 
late-time behavior and the number of periods of accelerated expansion. The 
concept of a quasi fixed point introduced in section 2.4 proves very useful for 
this purpose as it captures the asymptotic behavior of run-away solutions. 
This allowed us to calculate numerical values for cosmological observables at 
these (quasi) fixed points (given in tablelHI). We also discovered a new class of 
open universe solutions occurring in hyperbolic flux compactifications of M- 
theory. These undergo two periods of accelerated expansion where the first 
is transient while the second is eternal and asymptotes to zero. The quasi 
fixed point governing these solutions resides very close to the hypersurface 
of flat universes. Therefore all solutions approaching it are rather fiat at late 
times. 

The main motivation of our investigation is to describe the current uni- 
verse by interpreting the scalar field as a quintessence. As already noted by 
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many authors before, this model has the necessary quahtative feature that it 
generically induces accelerated expansion. The acceleration can be transient, 
eternal or eternally cyclic switching between acceleration and deceleration. 

Comparing the (quasi) fixed point solutions of our model with cosmo- 
logical measurements, however, needs some caution. Since the accelerated 
expansion is a relatively recent phenomenon in the history of the universe, 
it is clear that we cannot be arbitrarily close to a (quasi) fixed point yet. 
We expect that even in the case where our universe approaches such a point, 
the cosmological parameters computed in section 4 will receive corrections. 
Therefore it is hard to draw conclusions by comparing these quantities with 
observable data, as the (quasi) fixed points might only rule the universe some 
time in the future. 

A crucial ingredient missing in our setup is the component of (dark and 
ordinary) matter. Ref. [SH ES] list some options how this sector could be 
obtained from string theory, but it still remains a big challenge to derive the 
appropriate dark energy and matter from a fundamental theory simultane- 
ously. The inclusion of matter probably also modifies some of the quantities 
calculated at the (quasi) fixed points, although the qualitative features should 
remain the same. 

As investigated in ref. [34J, models coming from hyperbolic flux com- 
pactifications of M-theory with matter added by hand allow fine-tuned so- 
lutions, which are (marginally) compatible with cosmological data. Let us 
point out, that the double exponential potential derived from hyperbolic flux 
compactifications satisfies the tracker condition T = > 1 IHEI- For the 

potential (|7j) with /ci = — 1, F > 1 is always satisfied. It deviates significantly 
from unity as long as both exponential terms are of comparable size and ap- 
proaches r = 1 from above, once one of the exponential terms dominates. 
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This raises the hope that this model could give rise to a tracker behavior, 
which means that for a wide range of initial conditions the quintessence en- 
ergy density approaches and follows the background (radiation or matter) 
density before overtaking it and inducing a period of accelerated expansion. 

Another crucial point in the models with a higher- dimensional origin is 
the problem of a spontaneous decompactification of the extra dimensions. 
For the M-theory example of section 4, taking the age of the universe to be 
approximately 1.5 x 10^° years, and assuming that the solution describes the 
whole evolution except the big bang and primordial inflation, we find that 
the radius of the extra dimension should have increased by about a factor 
10'^. The current observation of three large dimensions is then compatible 
with a start value for r^l well below the Planck scale, even if one uses the 
naive estimate ttikk ~ '"int ^r the masses of Kaluza-Klein states. In the case 
of a hyperbolic internal manifold this should be a safe assumption, as these 
masses seem to be much higher than suggested by this naive estimate 

Our system also allows for an interpretation in terms of the very early 

universe. In this case matter effects are usually neglected. To realize a 

power law inflation scenario and solve the horizon and flatness problems 

without fine-tuning, one needs an attractive (quasi) fixed point with A > 1 

on the k = surface fZj . Comparing to table El we find that there are 

such points, viz. QFP2 and QFP7 for c > 3 and QFP5 for < c < 1. 

Solutions attracted to them are driven to = 1, even when k = ±1. The de 

Sitter fixed point, FP3, allows exponential expansion, but it is a saddle point 

and therefore the solutions need a drastic fine-tuning to produce inflation. 

Unfortunately none of those fixed points occur for scalar potentials obtainable 

from compactifications on maximally symmetric spaces with flux. 

""^""^In fact, for the work of ref. |17| which we generahzc, this was the main motivation. 
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While ten- and eleven- dimensional supergravities satisfy the strong en- 
ergy condition, p + 3p > 0, which forbids cosmic acceleration, an action 
of the form (jHl) only needs to satisfy the null energy condition, p + p > 0, 
which leaves some room for acceleration. However, if a theory is obtained 
by dimensional reduction of a theory obeying the strong energy condition, 
the lower-dimensional theory cannot have de Sitter vacua [SZI. This could be 
misinterpreted as a no-go theorem for accelerated expansion, but, as pointed 
out in ref. [HI, acceleration is nevertheless possible, if the internal mani- 
fold is time-dependent. However, in all known cases the acceleration is either 
transient, or eternal but going to zero asymptotically. In our tables (for com- 
pactifications, 1 < c < -\/3) this is reflected by the values of w for (quasi) fixed 
point solutions, which always satisfy and often saturate the bound w > 
set by the strong energy condition. This indicates how the strong energy 
condition satisfied by the higher-dimensional theory restricts the amount of 
acceleration obtainable in time-dependent compactifications, by requiring it 
to converge to zero for late times. Note that this fits nicely with the fact 
that the internal space always decompactifies in this limit. Therefore it 
is natural that the asymptotic behavior can be understood in terms of the 
higher-dimensional theory. 

Let us finally discuss the main approaches for avoiding the consequences 
of the strong energy condition valid in the higher- dimensional theory and 
for stabilizing moduli. Constructions involving non-perturbative effects of 
string/M-theory, like branes, orientifolds, and instantons can produce infla- 
tion and stabihze the moduli jSHj, but seem to require a significant discrete 
^^This observation agrees with the more general analysis [2Sj which argued that com- 
pactifications leading to an accelerated expansion of the 4-dimensional space-time either 
decompactify at late-times or end in a big crunch singularity. 
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fine tuning. Therefore it remains important to investigate various otlier di- 
rections. For example, tfie dimensional reduction on group manifolds also 
leads to multi exponential potentials allowing accelerating cosmologies [22] • 
One promising strategy to avoid the strict limits on the amount of accelera- 
tion obtainable in compactifications is to consider generalized compactifica- 
tions, where one twists the theory by a symmetry of the equations of motion 
which is not a symmetry of the action. This leads to lower-dimensional the- 
ories which do not have an action j3U].^^ It would be interesting to explore 
whether such constructions can lead to double exponential potentials with c 
being outside the range 1 < c < -\/3 valid for maximally symmetric spaces. 
Another interesting direction was pointed out by the authors of refs. |42j . 
who argue that the 'central region' in the moduli space, which cannot be 
mapped to a perturbative regime by duality transformations, might give vi- 
able cosmologies without fine-tuning. One particular approach for exploring 
this region is the inclusion of states which become light at special points in 
the moduli space The effects of such extra states improve on the moduli 
problem, seem to be important for vacuum selection, and help to get cosmic 
acceleration H^l • 
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